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Abstract. The hydrodynamic limit of the kinetic Cucker-Smale flocking model 
H is investigated. The starting point is the model considered in [19], which in 

C ' 1 addition to the free-transport of individuals and the Cucker-Smale alignment 

£SJ operator, includes a strong local alignment term. This term was derived in [20] 

as the singular limit of an alignment operator due to Motsch and Tadmor [25] . 
The model is enhanced with the addition of noise and a confinement potential. 
The objective of this work is the rigorous investigation of the singular limit 
corresponding to strong noise and strong local alignment. The proof relies 
on the relative entropy method and entropy inequalities which yield the ap- 
propriate convergence results. The resulting limiting system is an Euler-type 
flocking system. 
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1. Introduction 

Mathematical models aimed at capturing parts of the nocking behavior exhib- 
ited by animals such as birds, fish, or insects, are currently receiving widespread 
I attention in the mathematical community. Many of these models have sprung out 

in the wake of the seminal paper by Cucker and Smale [9] . The typical approach 
. is based on particle models where each individual follows a simple set of rules. To 

date, the majority of studies on flocking models have been on the behavior of the 
particle model or the corresponding kinetic equation. 

For practical purposes, if the number of individuals in the flock is very high, it 
might be desirable to identify regimes where the complexity of the model may be 
reduced. This paper is a modest contribution to this complex task. The starting 
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point for our study is the following kinetic Cucker-Smale equation on E d x M. d x (0, T) 

ft +v ■ VJ + div v (/£[/]) - V,* • V„/ - <rA v f + /3 div„(/(« - u)). (1.1) 

Here, / := f(t,x,v) is the scalar density of individuals, d > 1 is the spatial dimen- 
sion, (3, a > are some constants and $ is a given confinement potential. The 
alignment operator L is the usual Cucker-Smale (CS) operator, which has the form 

L[f}= [ I K(x,y)f(y,w)(w-v)dwdy, (1.2) 

JR d JR d 

with K being a smooth symmetric kernel. The last term in (1.1) describes strong 
local alignment interactions, where u denotes the average local velocity, defined by 

J Rd f dv 

This strong alignment term was introduced in [20] as the following singular limit 
of the Motsch-Tadmor (MT) alignment operator 

hmL / = lim ^* r rj~ t r^— ~, = {u- v). (1.3) 

] Rd J Rd (j)(x-y)f{y,w)dwdy 

Since the MT term is relatively new in the literature, a remark on its purpose seems 
appropriate. The MT operator was introduced in [25] to correct a deficiency in the 
standard Cucker-Smale model. Specifically, since the CS operator L[] is weighted 
by the density, the effect of the term is almost zero in sparsely populated regions. 
The MT operator instead weights by a local average density. To arrive at (1.1), 
the rationale is to let the MT operator govern alignment at small scales and the 
CS operator the large scales. Our equation (1.1) is then obtained in the local limit 
(1.3) which seems appropriate at the mesoscopic level. 

Since the kinetic equation (1.1) is posed in 2d+l dimensions, obtaining a numer- 
ical solution of (1.1) is very costly. In fact, the most feasible approach seems to be 
Monte-Carlo methods using solutions of the underlying particle model with a large 
number of particles and realizations. Consequently, it is of great interest to deter- 
mine parameter regimes where the model may be reduced in complexity. The goal 
of this paper is to study the singular limit of (1.1) corresponding to strong noise 
and strong local alignment, that is a, — > oo. More precisely, we are concerned 
with the limit e — > in the following equation: 

n+v-v x r+div v (rL{r})~v x <s>-v v r = ~A v r+ 1 dw v (r(v-u')). (1.4) 

e e 

This scaling can alternatively be obtained from (1.1) by the change of variables 

x = ex, t — et, 

and assuming that K{x, y) = eK(x — y). 

In the remaining parts of this paper, we shall establish with rigorous arguments 
that 

> "-^(*.°oi 2 
f ^ g(t,x)e 

where g and u are the e — > limits of 

g £ = [ p dv, g £ u e = [ f e v dv. 

As a consequence, we will conclude that the dynamics of / is totally described by 
the following Euler-Flocking system 

Qt + dw x (gu) = 0, (1.5) 
(gu)t + div x (gu®u) + S7 x g= K \x , y) g(x) g(y)[u(y) - u(x)\ dy - gV x <& . (1.6) 
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The result will be precisely stated in Theorem 3.1 with the proof coming up in 
Section 4. The proof is established via a relative entropy argument providing in 
addition a rate of convergence in e. The relative entropy method relies on the "weak- 
strong" uniqueness principle established by Dafermos for systems of conservation 
laws admitting convex entropy functional [11] (see also [12]). It has been successfully 
used to study hydrodynamic limits of particle systems [16, 21, 24, 27]. In our case, 
the result will be somewhat restricted as we need the existence of smooth solutions 
to (1.5) which we only know locally in time. 

Remark 1.1. An important issue in the study of Cucker-Smale type equations is 
whether a given model leads to flocking behavior. By flocking, it is usually meant 
that the velocity u(x, t) converges, for large t, to a constant velocity u. Here, we 
note that flocking can only occur in (1.5)-(1.G) if the confinement potential and 
the pressure are in balance. Indeed, if (g,u) (with u constant) solves (1.6), then 
\7g = — gV$ and hence g — j e K 'i da . e~ j '. We thus see that flocking only occurs 
when the density is very diluted and, in particular, does not have compact support. 

1.1. Formal derivation of (1.5) - (1.6). For the convenience of the reader, let 
us now give the formal arguments for why (1.5) - (1.6) can be expected in the limit. 
First, we note that to obtain an interesting limit as e — > in (1.4), the right-hand 
side should converge to zero 

A v f + div v (f e (v-u e ))^0. 

If this is the case, the limit / can only have the following form 

_ \v-U,{t,3s)\ 2 

f -> f{t,x,v) = g(t,x)e 

Hence, it seems plausible that the evolution of / (in the limit) can be governed by 
equations for the macroscopic quantities g and u alone. 

To derive equations for g and u, let us first integrate (1.4) with respect to v 

g\ + div(/u e ) = 0. (1.7) 

Hence, by assuming the appropriate converge properties and passing to the limit, 
we obtain the continuity equation (1.5). 

To formally derive (1.6), let us multiply (1.4) by v and integrating with respect 
to v to obtain 

(g e u e ) t + div x ( / {v®v)f t dv 

(1.8) 

K(x, y)g e {x)g'(y)(u e (x) - u e {y)) dy + e V*$ = 0. 

Passing to the limit in (1.7) and (1.8) (assuming that f e — > ge 2 ( g e — > g and 
u e — > u), we get: 

(gu) t + div x [ g j (v ® v)e~~ ^~ dv 

(1.9) 

K(x, y)g(x)g(y)(u(x) - u(y)) dy + gV x $ = 0. 
By adding and subtracting u, we discover that 

. 1 



(v <g) v)e 2 dv = (u ® u)e 2 + (v — u) <£> (v — u)e 2 dv 
= u ® u + 1. 

Inserting this expression in (1.9) gives (1.6). □ 



4 



KARPER, MELLET, AND TRIVISA 



Organization of the paper: The rest of this paper is organized as follows: In 
Section 2, we recall some existence results for the kinetic flocking model (1.1) (these 
results were proved in [I 1 )]) and for the Euler-flocking model (1.5)-(1.6) (a proof of 
this result is provided in Appendix A). In Section 3 we present our main result, 
which establishes the convergence of weak solutions of the kinetic equation (1.4) to 
the strong solution of the Euler-flocking system (1.5)-(1.6). The proof of the main 
theorem is then developed in Section 4. 

2. Existence theory 

The purpose of this section is to state some existence results upon which our 
result relies. More precisely, the proof of our main result (convergence of (1.4) to 
(1.5)-(1.6)) makes use of relative entropy arguments which require the existence 
of weak solutions of (1.4) satisfying an appropriate entropy inequality, and the 
existence of strong solutions to the Euler-Flocking system (1.5)-(1.6) satisfying an 
entropy equality. Note that the latter result will be obtained only for short time. 

Since entropies play a crucial role throughout the paper, we first need to present 
the entropy equalities and inequalities satisfied by smooth solutions of (1.4) and 
(1.5)-(1.6). 

2.1. Entropy inequalities. Solutions of (1.1) satisfy an important entropy equal- 
ity, which was derived in [!!)]: We define the entropy 

12 



Hf)= / flogf + f^ + fQdvdx (2.1) 



and the dissipations 

Dx(f)= I -AV v f-f(u-v)\ 2 dvdx, 

2d (2 2) 

AaCf) = x / / / / K(x,y)f(x,v)f(y,w)\v-w\ 2 dwdydvdx. 

Z JS. d JR d JR d JR d 

The latter is the dissipation associated with the CS operator L[-}. There holds: 

Proposition 2.1. Assume that L is the alignment operator given by (1.2) with K 
symmetric and bounded. If f is a smooth solution of (l-4)> then f satisfies 

d t F{f)+ l D 1 {f) + D 2 (f) 
e 

K(x, y)f(x, v)f(y, w) dw dy dv dx, (2.3) 

with J-'(-), -Di(-), -D2G) given by (2.1) and (2.2). Furthermore, if the confinement 
potential $ is non-negative and satisfies 

e-*< s > dx < +00, (2.4) 
then there exists C, depending only on \\K\lao, <3? and J fo(x,v) dxdv, such that 

+ \ f [ K(x,y)g(x)g(y)\u(x)-u(y)\ 2 dy dx < CeT(f(t)). 

1 JR d JM d 

The first inequality (2.3) shows that the nonlocal alignment term is responsible 
for some creation of entropy. The second inequality (2.5) shows that this term can 
be controlled by D\{f) and the entropy itself. This last inequality will play a key 
role in this paper. 
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The Euler system of equations (1.5)-(1.6) also satisfies a classical entropy equal- 
ity. More precisely, if we define 

$(q,u) = / g— dx + / g\og g + g&dx, 
Jm d 2 J R d 

then any smooth solution of (1.5)-(1.6) satisfies 

d t £(Q,u) + -- I I K(x,y)g(x)g(y)\u(x) - u{y)\ 2 dydx = 0. 
1 JR d Jm d 

Note that the entropy & and T are related to each other by the relation 

Furthermore, we have the following classical minimization principle (consequence 
of Jensen inequality): 

#{Q,u)<T{f), iff?= / fdv, gu = J vfdv. (2.6) 

This relation will be important in the upcoming analysis when considering the 
relative entropy of solutions to the kinetic equation (1.4) and solutions to (1.5) - 
(1.6). 

2.2. Global weak solutions of the kinetic equation. The existence of a weak 
solution for (1.1) is far from trivial because of the singularity in the definition of u. 
We will say that a function / satisfying 

/ G C(0,T;i 1 (M 2d ))n J L oo ((0,T) x R 2d ), (\v\ 2 + $(*))/ G L°°(0, oo; L 1 (R 2d )), 

is a weak solution of (1.1) if the following holds: 

/ -fit - ^/V.V + /VxSV^ - fL[f\V v 1> dvdxdt 
+ / <TV v fV v ip - f3f(u - v)V v tp dvdxdt (2.7) 
f°i/)(0, •) dvdx, 



for any -0 G C^°([0,T) x E 2d ), where u is such that j = gu. 

Remark 2.2. Note that the definition of u is ambiguous if g vanishes (vacuum). We 
resolve this by defining u pointwise as follows 

r jjpt) xt \^ Q 

u{x,t) = I g{x,t) ^ ' Jr- (2.8) 
{ if g(x,t) = 

This gives a consistent definition of u as can be seen from the bound 

3< (J \v\ 2 f(x,v,t)dv^j 1 g 1 ' 2 , 

yielding j = whenever g — and so (2.8) implies j — gu. 

Note also that u does not belong to any LP space. However, we have 



\uf\ dxdv < ||/||Loo (R 2d) / \v\ f(x,v,t)dvdx, 

JR 2d 

so the term uf in the weak formulation (2.7) makes sense as a function in L 2 . 

The existence result we shall utilize in this paper was obtained as the main result 
in [19] and is recalled in the following theorem: 
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Theorem 2.3. Assume that L is the alignment operator (CS) given by (1.2) with 
K symmetric and bounded. Assume furthermore that fo satisfies 

f eL oo (R 2d )nL 1 (M 2d ), and {\v\ 2 + $(z))/ € L 1 (M 2d ). 

Then, for all e > 0, there exist a weak solution f e of (1-4) satisfying 

Hf c {t)) + [ -Di(f e ) + W) ds < J-(/o) + Ct, (2.9) 
Jo e 

where the constant C depends only on \\KWoo, $ and J fo(x,v) dx dv. Furthermore, 
if satisfies (2.4), then f e also satisfies 

Hf e («)) + ~ f D^Dds+l [ f K(x,y)g*{x)g\y)\u*{x)-u*{y)\ 2 dydxd S 

JO 1 Jo jR2d 

<T(f )+Ce f F{r{s))ds (2.10) 
Jo 

for all t > 0. 

2.3. Existence of solutions to the Euler-Flocking system. As usual with 
relative entropy methods, our main result will state that the solutions of (1.4) 
converge to a strong solution of the asymptotic system (1.5) - (1.6), provided such 
a solution exists. It is thus important to prove that (1.5) - (1.6) has a strong 
solution, at least for short time. This is the object of the next theorem (which we 
state in the case d = 3): 

Theorem 2.4. Let (po,u ) € H S (R 3 ) with s > 5/2 and p (x) > in R 3 and 
assume that V x $ € H S (M. 3 ). Then, there exist T* > and functions (g,u) £ 
C{[0,T*];H s (R 3 )) n C 1 ((0,T*);H s - 1 (R 3 )), p(x,t) > 0, such that (g,u) is the 
unique strong solution of (1.5) - (1.6) fort € (0, T*). Moreover, (g,u) satisfies 
the equality 

d t S{g,u) + \( [ K(x,y)g(x)g(y)\u(x)~u(y)\ 2 dydx^O. (2.11) 

1 Jm d Jm d 

Since the proof of this theorem is rather long and independent of the rest of the 
paper, we postpone it to the appendix. 

Note in particular that the condition s > 5/2 implies that the solution satisfies 

ue L°°{[0,T*};W hoo (R 3 )). (2.12) 
Furthermore, dividing the momentum equation by /?, we also get: 

V^logp € L°°([0,T*] x R 3 ). (2.13) 
These two estimates is all the regularity we will need in our main theorem below. 

3. Main result 

With the existence results of the previous section, we are ready to state our main 
result concerning the convergence of weak solutions of (1.4) to the strong solution 
(g, u) of the Euler-flocking system (1.6)-(1.5) as e — > 0. 

Theorem 3.1. Assume that: 
(1) fo is of the form 

_ g Q (x) _ i"o(*)-H 2 
J° (2^/2) 6 

with 

fo G L°°(M. 2d ) fl L 1 (R 2d ), and (\v\ 2 + *(x))/ G L 1 ^). 
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(2) / e is a weak solution of (1.4) satisfying the entropy inequality (2.10) and 
with initial condition / c (0, •) = fo(-). 

(3) T* > is the maximal time for which there exists a strong solution (g, u) 
to the Euler system of equations (1.5) - (1-6), with g$ = J Rd fo dv and 
qqUq = J Rd fav dv and satisfying (2.12) and (2.13) (Theorem 2.4 gives in 
particular T* > if po and Uq are regular enough). 

There exists a constant C > depending on 

F(fo), ||-K1U°°, T *, IMU~(0/r*;Wi.«>(R<*))! and ll vl °gHlL»((o,r«)xr)> 
such that 

— \u e - u\ 2 + I" ^— ^ dz dxdt (3.1) 
2 ./„ z 



T 

JR d 



K(x, y)g e (x)g e (y) [{u e (x) - u(x)) - (u e (y) - u(y))} 2 dxdydt 



1 

+ 2 
<Cyfe, 
where 

g e = [ f e dv, g e u e = [ vp dv. 

Moreover, any sequence of functions satisfying (3.1) satisfies: 

f e ^ ge'^ 1 ^ a.e and Lj oc (0, T* ; l}(m. d x R d )), 

g e ^> g a.e and Lj oc (0,T*;Z ( 1 (R <i )), 

g e u e ^ gu a.e and Lj oc (0, T*;L 1 (R d )), 

g e \u e \ 2 ^ gu 2 a.e and Lj oc (0, T*;L 1 (R d )). 

This theorem will be a direct consequence of Proposition 4.1 below, the proof of 
which is the object of Section 4. 

4. Proof of Theorem 3.1 

To reduce the amount of notations needed in the proof of Theorem 3.1 it will be 
preferable to write the Euler-Flocking system in terms of the conservative quanti- 
ties. In our case, the conservative quantities are the density g and the momentum 
P = gu. If we denote 



u = 



p 



we can rewrite the system (1.5)-(1.6) as 

U t +dW x A(U)=F(U). (4.1) 
The flux and source term are then given by 

A (C0 =(i °), F(V) = ( eP _J_^)' 

where we have introduced the notation g = J Rd K(x, y)g{y) dy. The entropy E 
corresponding to (4.1) reads 

p2 

E(U) = — + elogQ + Q$, 
2g 

and the relative entropy is the quantity 

£{V\U) = E(V) - E{U) - dE(U){V - U), 
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where d stands for the derivation with respect to the variables (g,P). 
For the system (1.5)-(1.6), a simple computation yields 

-dE(U)(V-U) = -^ +l J + 1 + (*Z° p 

<t 1 7 / 1 2 nlni* 

+ (q — ?)(log|0+ 1 + $) + Qu - quv. 



2 2 

Conseqently, the relative entropy can alternatively be written 
£{V\U) = E(V) - E{U) - dE(U){V - U) 

Id 2 Id 2 
= q~Y~ ~ Q-^ + Q^ogq- glogg + ^(q- g) 

q\u\ 2 o\u\ 2 , .., 9 

+ ~ ^y- + (e- g)0oge + i + $) + ^ 2 - ?™ 

|2 



(4.2) 



It! — u\ 

= g J 2 +p(glg), 

where we have introduced the relative pressure 

p(?|g) = glogg-£log£+(g-g)(log£> + l) = / dz. 

J g z 

Note that the relative pressure controls the L 2 norm of the difference 

p(q\g) > ^-mml -}-,-}-} (q(x) - g(x)) 2 . (4.3) 



2 [ q(x) ' q(x) 

With the newly introduced notation, Theorem 3.1 can be recast as a direct 
consequence of the following proposition: 

Proposition 4.1. Under the assumptions of Theorem 3.1, let 

U=( Q 
\gu 

denote the strong solution to the Euler system of equations (1.5)-(1.6) and let 
U € = (?.) , Q e = ! r dv, g e u* = [ f% dv, 

\P U J JR d JR<> 

be the macroscopic quantities corresponding to the weak solution of the kinetic equa- 
tion (1.4). 

The following inequality holds: 

£{U £ \U){t) dx 



1 



2 Jo 

< C 



K{x,y)g e {x)g t {y) [(u e {x) - u(x)) - (u e (y) - u(y))] dxdyds 
f f £{U e \U)dxds + C«Ji. 

JO JR d 



(4.4) 



The proof of this proposition relies on several auxiliary results which will be 
stated and proved throughout this section. At the end of the section, in Section 4.6 
we close the arguments and conclude the proof. However, before we continue, let 
us first convince the reader that Proposition 4.1 actually yields Theorem 3.1. 
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Proof of Theorem 3.1. Let us for the moment take Proposition 4.1 for granted. 
Then, the main inequality (3.1) follows from (4.5) and Gronwall's lemma. 

We now need to show that (3.1) implies the stated convergence. First, we note 
that the entropy estimate (2.9) implies that f e is bounded in LlogL and thus 
converges weakly to some /. 

Next, in view of (4.2) and (4.3), the main inequality (3.1) yields 

[ [ mini — , -1 \p £ - p\ 2 dxdt — >0 
Jo Jm« \p e P) 

[ [ p e \u e -u\ 2 dxdt — >0. 
Jo Jm d 

In particular, we get: 
/ \p e — p\ dx = [ mini— , — I max{p e , p} 1 ^ 2 \p t — p\ dx 

Jr* ' 1 Jr* Ip £ p) 



and 



1 1 



< ( / min <{ - ^ \p e — p\ 2 dx ) ( / ma,x{ p e , p}dx 



< / min 



1 1 

P e ' P 



\ 1/2 

\p e -p\ 2 dx) (2M) 1 / 2 



Hence, we can conclude that 

g e ^ q a.e and L\ oc {0, T* ; L 1 (M. d )) , 
Similary, we see that 

\g e u e — gu\ dx < / \g e (u e — u)\ + \(g e — g)u\ dx 



<M5 g e \u e -u\ 2 dx 



1 1 

P 6 ' P 



P - P\ 



v 1/2 / „ s 1/2 

|2 dx J \J max{p e , p}u 2 dx J 



Consequently, also 



^ gu a.e and h\ oc (0, T* ; L 1 (R d )) . 



g u 

Moreover, by writing 

p e u e2 - pu 2 = p e (u e - u) 2 + 2u(p e u e - pu) + u 2 (p- p c ) 
we easily deduce 



g 



VI 2 ^ Qu 2 a.e and L} oc (0,T*;L 1 (R d )). 



At this stage, it only remains to prove that / has the stated maxwellian form. 
For this purpose, we first send e — > in the entropy inequality (2.9) and use the 
convergence of g e u e and p e |u e | 2 to obtain 

lim / f^-+f \ g f dx+l f f K(x,y)g(x)g(y)(u(y)-u(x)) 2 dxdydt 

e ^"J«.d & & Jo JR d 

< £(qo,u ), 

(4.5) 
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-(v-uq) 2 

where we have used that /o = ^ 2 ^° d/2 e = to conclude the last inequality. Next, 
we subtract the entropy equality (2.11) from (4.5) to discover 

r v 2 r u 2 

o<iim/ r- +r log r dxdv~ B +e \ oged x<o 1 (4.6) 

where the first inequality is (2.6). By convexity of the entropy, we conclude that 

(5hr)l 

which concludes the proof of Theorem 3.1. □ 

4.1. The relative entropy inequality. The fundamental ingredient in the proof 
of Proposition 4.1 is a relative entropy inequality for the system (1.5) - (1.6) which 
we will derive in this subsection. However, before we embark on the derivation of 
this inequality, we will need some additional identities and simplifications. 

First, we recall that E is an entropy due to the existence an entropy flux function 
Q such that 

djQi(U) = Y J d j A ki {U)d k E{U), i,j = l,...,2. (4.7) 

k 

We then have 

E(U) t +div x Q{U) = (-QP + QP)-. (4.8) 

Q 

Since the confinement potential term 3>q in the entropy is linear, it does not 
play any role in the relative entropy. It is thus convenient to introduce the reduced 
entropy functional 

E{U) = — + Q\0gQ. 

2q 

This allows us to treat the contribution of <&q as a forcing term (which is part of 
the F(U)) in the proof of Proposition 4.2. The reduced entropy flux Q(U) is then 
defined by 

djQiiU) = Y,djA ki {U)d k E{U), i,j = l,...,2. (4.9) 

k 

We note that Q satisfies 

div x Q(u) = dE(U) (div x A{U)) , (4.10) 
while the total entropy flux, defined by Q(U) — Q(U) + P$, satisfies 

div x Q(U) = (div x A(U) + QvJ) dE(U). 

We shall also need the relative flux: 

A(V\U) = A(V) - A(U) - dA(U)(V - U), 
where the last term is to be understood as 

[dA(U)(V - U% = dAi(U) >{y-U), i = l....,d. 
The key relative entropy inequality is given by the following proposition: 

Proposition 4.2. Let U = ( be a strong solution of (4.1) satisfying (2.11) 

and let V = ^ ^ be an arbitrary smooth function. The following inequality 
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holds 
d 

dt 



£(V\U) dx 



< 



K(x, y)q(x)q(y) [(v(x) - u(x)) - (v(y) - u{y))] dxdy 



d t E(V)+QV x <P + ^ I K{x,y)q{x)q{y)[v{x)-v(y)] 2 dy 



dx 



V x {dE{U)) : A{V\U) dx 
dE(U)[V t + divA(V)-F(V)] dx 

K(x, y)q{x)(g(y) - q(y))[u(y) - u(x)][v(x) - u(x)] dxdy 



Remark 4.3. When K = 0, such an inequality was established by Dafermos [11] for 
general system of hyperbolic conservation laws. 

In order to prove Proposition 4.2, we will need the following lemma (see Dafermos 
[ ])• 

Lemma 4.4. The following integration by parts formula holds 
d 2 E{U) (di V:r A{U)) {V - U) dx 



(dA(U){V -U)) : (v x dE(U)} dx, 



(4.11) 



where : is the scalar matrix product. 

Proof of Lemma 4-4- The proof of this equality can be found at several places in the 
literature. For the sake of completeness we recall its derivation here (we follow [ ] 
[p. 1812]). 

Differentiating (4.9) with respect to Ui, we obtain the identity 

d d 

Y / d i d kE(U)d J A kt {U) = didjQi(U) -J2 d i d 3 A ki(U)d k E{U). 



k=l 



k=l 



Using this identity, we calculate 
d 2 E(U) (div* A{U)) (V - U) 



d I d \ 
lij=l \k=l I 



dxi 



(Vi - Ui) 



3U - BIT- 
^didjQiffl-^lVi-U,)- didjAkiWdkEW-^m-Ut) 

£ diy(diQ(U)(Vi - U x )) - div - i d iMU)(Vi - Ui)) d k E(U) 



kl 



+ J2(V x V t -V x Ui) 



J2 d iA k (U)d k E(U) ~ diQ(U) 
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Now, we observe that (4.9) implies that the last term is zero. Thus, we can conclude 
d 2 E{U) (div x A(U))(V-U) 

= div x (dQ(U)(V - U)\ - div x (dA(U)(V - U)) dE(U) 

Integrating this identity over M. d yields (4.11). □ 
Proof of Proposition 4-2. First of all, we recall that 

£(V\U) = E(V) - E(U) - dE{U){V - U) 
= E(V) - E{U) - dE{U){V - U). 

We deduce: 

S(V\U)dx= [ d t E(V)-dE{U)U t -d 2 E(U)U t (V -U) 
"* Jm. d Js. d 

-dE(U)(V t -U t ) dx 
d t E{V) - d 2 E(U)U t (V — U) — dE(U)V t dx 
:=h+I 2 +h. 
Since d 2 E(U) = d 2 E(U), formula (4.11) provides 

d 2 E(U) [div x A{U) - F(U)] {V - U) dx 



(4.12) 



V x dE(U) : dA{U){V - U) dx - / d 2 E(U)F(U)(V - U) dx. 
By adding and subtracting, and integrating by parts, we find 
I 3 = - [ dE(U) [V t + div A(V) ~ F(V)] dx 

- [ (v x dE(U)) : A(V) + dE(U)F(V) dx 

jR d ^ 1 

Consequently, 

h + h = ~ [ V x {dE{U)) : A{V\U) + dE(U) [V t + divA(V) - F(V)] dx 

V x (dE{U)) : A{U) dx (4.13) 

- / d 2 E{U)F{U){V -U) + dE{U)F{V) dx := J x + J 2 + J 3 . 

Now, using (4.10), we get: 

J 2 = - / V x (dE(U)) : A(U) dx 

/ Rd , (4-14) 

= / dE{U) d\v x A{U) dx= div x Q{U) dx = 0. 

JM d JR d 

It only remains to compute the term J3, which is the only non standard term, 
since it includes all the contributions of the forcing term F(U). We now insert our 
specific expression of E and U. A simple computation yields: 

dm=( d A { !!)) = (-^ +1 ° sg+1 ), <?e(u)={\ tV 



d P E(U) 



ij 
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Let us also introduce two functions (q, Q) such that V = [q, Q] T and define v = ®. 



-J 3 = / d 2 E(U)F(U)(V -U) + dE{U)F(V) dx 

jR d 

= [ -^(qP-qP- g\7 x <S>) (q-g) + - (gP - qP - g\7 x <S>) (Q - P) 
J R d g z \ ) g\ ) 

+ ^ {qQ -qQ- ?V x $) dx 

= [ - {qP- QP) (-- P ) + P (qQ - qQ) - gV x $ dx 

= / q (qu - B u ) {v — u) +u[qQ — qQ) - QV X <& dx 

jR d v ' 

which yields 

-■h = I / K(x,y)q(x)g(y)[u(y) -u(x)][v(x) -u(x)] 

JR d JR d 

+ K(x, y)q(y)q(x)u(x)[v(y) - v(x)] dydx. 
= / K(x,y)q(x)q(y)[u(y) -u(x)][v(x) -u(x)] dxdy 

+ J J K{x,y)q{x)q{y)[v(y) - v(x)][u(x) - v(x)} dxdy ^ ^ 

+ / K(x,y)q(x)q(y)[v(y) -v(x)]v(x) dxdy 

JR d JM d 

+ / K(x,y)q(x)(g(y)-q(y))[u(y)-u(x)][v(x)-u(x)]dxdy 

JR d JR d 

qvV x & dx. 



/R d 

Using the symmetry of K, we see that the first two terms can be written 
/ / K {x,y)q(x)q{y)[u(y) - u(x)][v(x) - u(x)] dxdy 

JR d JR d 

+ / K{x,y)q{x)q{y)[v{y) - v{x)] [u(x) - v(x)] dxdy (4.16) 

JR d JR d 



R d JR d 

K(x,y)q(x)q(y)[u(y) - u(x)][(v(x) - u(x)) - (v(y) - u(y))] dxdy 



K(x,y)q(x)q(y)[v(y) - v(x)][(u(x) - v(x)) - (u(y) - v(y))] dxdy 

= nl I K ( x 'y)q( x )q(y)i( v ( x ) - u ( x )) - ( v (y) - u (y))} 2 dxd v- 

L JR d JR d 

From the symmetry of K, we also easily deduce 

K(x,y)q(x)q(y)[v(y) - v(x)]v(x) dxdy 

1 f f (4 - 1? 

= -- / / K(x,y)q(x)q(y)[v(x) - v(y)} 2 dxdy. 

z JR d JR d 
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Hence, by setting (4.17) and (4.16) in (4.15), we discover 



J* = 



' ' K(x,y)q(x)q(y)[v(x) - v(y)] 2 dxdy + I qvV x § dx 

K(x, y)q(x)q(y) [(v(x) - u(x)) - (v(y) - u{y))] 2 dxdy 
+ I K(x, y)q(x)(g(y) - q{y))[u{y) - u(x)][v(x) - u{x)] dxdy 



We conclude the proof by combining the previous identities (4.14), (4.13) and (4.12). 

□ 

In our proof Proposition 4.1, we will use the following immediate corollary of 
Proposition 4.2: 

Corollary 4.5. Let / e be a weak solution of (1.4) satisfying (2.10) and let 
U e = (e» c ,£»V), with g e = / fdv, g e u e = [ vfdv. 

JS. d JR d 

Let U = (g, gu) be the strong solution of (4.1) satisfying (2.11). Then the following 
inequality holds: 

| / E{W\U) dx 
at . U d 



< 



K(x, y)g e (x)g e (y) [(u e (x) - u(x)) - (u c (y) - u(y))} dxdy 
d t E(U e ) + P e V x $ + l 



dx 



K(x,y)g'(x)g*(y)[u*(x)-uC(y)] 2 dy 
V x {dE{U)) : A{W\U) dx 

dE{U)[U f t +div A{U e ) - F(U e )] dx (4.18) 
K(x,y)g e (x)(g{y) - g e (y))[u(y) - u(x)][u e (x) - u{x)\ dxdy. 



In order to deduce Proposition 4.1 from this corollary, it remains to show that 

(1) The first term in the right hand side in (4.18) is of order e when integrated 
with respect to t (Lemma 4.6). 

(2) The second term is controlled by the relative entropy itself (in fact we will 
show that the relative flux is controlled by the relative entropy, see Lemma 
4.7). 

(3) The third term is of order 0(e) (Lemma 4.8). 

(4) The last term can be controlled by the relative entropy (Lemma 4.9). 
The rest of this paper is devoted to the proof of these 4 points. 

4.2. (1) The first term. 

Lemma 4.6. Let f c be the weak solution of (1.4) given by Theorem 2.3 and let 



U £ = (g e ,g e u e ), with g e 



fdv, g e u e 



vf e dv. 



Then 
ft 



JR d 



d t E(U')+P e V x <f> 



K(x,y)g*(x)g i (y)[u*(x)-u*(y)} 2 dy 



dx ds 



<Ce f Hf(s))ds 
Jo 
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for all t > 0. 

Proof. First, we write 

/ d t E(U e ) ds 

JR d 

' E{U e )(t) - E(U Q ) dx 

[E{U%t)-P{f%t)] + [Hf e )(t)-f(f )] + [?(f )-E(U )]. 
The well-preparedness of the initial data gives 

/ F(fo) - E(U ) dx = 0, 

JR d 

and (2.6) implies 

f E{W){t)-T{f){t)dx<Q. 
Jm d 

Finally, we deduce (using (2.10)) 

r d t E(U £ )ds< I F(f){t)-JF{f )dx 

<f Hf)(t)- Hfo) dx- f p € (t)$- po^dx 

JR d JWL d 
1 '* 



JM d 



<-~ I / / K(x,y)g e {x)Q e {y)\u e (y)-u c (x)\ 2 dydxds 

1 JO JM d JR d 

- [ [ P'VQdxds + Ce f F(f e (s))ds 

Jo JR d JO 



□ 

4.3. (2) Control of the relative flux. We note that \ V x dE(U) \ is bounded in L°° 
by |M|l°o( ,t* ; w 1 >°°), ||Vlogp||ioo, so the second term in (4.18) will be controlled 
if we prove the following standard lemma: 

Lemma 4.7. The following inequality holds for all U , V : 

[ \A{V\U)\dx < [ £{V\U) dx 

JR d Jw d 

Proof. A straightforward computation gives 
dA(U)(V-U) 



(g-g) 



Q-P \ 

y P®P+]P®{Q-P) + \(Q-P)®P+(q- Q)) 1 

and using the fact that P = gu and Q = qv, we get: 

-Q <8 Q - -P <8> P - -P <8 (Q - P) - -(Q - P) <8> P + - ~ - P <8 P 
q g g g g z 

— qv ® v + gu ® u — u ® qv — qv ® u + (q — g)u ® u 
= g(i> — u) ® (v — u). 

Since all the other terms in -A(-|-) are linear, we deduce 

A ( V W = (q(v-u)®(v-u) y 
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This implies 

/ |-A(F|Z7)I dx= [ q\v-u\ 2 dx< I £{V\U) dx, 
which concludes our proof. 

4.4. (3) Kinetic approximation. 



□ 



Lemma 4.8. Let U be a smooth function, let f e be a weak solution of (1.4) satis- 
fying (2.10) and define 



f e dv, g e u e = / vf dv 



U c = (g e , g e u e ), with g e 



There exists a constantC depending onT, ||'u||l oo (o,t*;W 1 >° o ) an d ||Vlogp||L°o such 
that 



dE{U) [W t + div A(W) - F{U e )] dx 



10 JR d 

Proof. By setting := <fi(t,x) in (2.7), we see that 

g e t +div(g e u e ) = in D'([0,T) X ft). 

Setting <f) :— v^(t,x), where "J is a smooth vector field, we find that 
(£>V) t + dW x (g e U e (g> u e ) + \7 x g e 

K{x,y)g t {t,x)g t {t,y){u t {x) - u e (y)) dy + g e V x $ 



(4.19) 



(4.20) 



= div^ J (u e <g u e - v <g v + I) f dv, 
in the sense of distributions on [0,T) x ft. Hence, we have that 
f f dE(U) [U e t + div A(U e ) -F(U e )} dxdt 

JO JR d 



< 



JR d 



\V x dE(U)\ 



(u e <g u € - v <g v + I) f e dv 



dxdt 



(4.21) 



< C 



JR d 



(u € <g u e - v <g> v + I) / c dv 



dxdt, 



where the constant C depends on | |w| \l°°(o,t*-,w 1 <°°) an d 1 1 V log To conclude, 

we have to prove that the righthand side can be controlled by the dissipation. As 
in [23], we calculate 

(u e (g) u e - v ® v + l)f e dv 

(w £ <g> (u e - v) + (u e - v) ® v + I) / £ dv 

u e y/T e ® ((u e -f) V / 7 7 -2V„ v /7 7 ) +« £ 8V„f 

+ (V ~v)y/f* -2V 1 , V // 7 ) ®v^ + V v f e ®v + lf e dv. 
Using integration by parts, we see that 



u e (g V„f dv = 0, 



V w / £ <g v dv 
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By applying this and the Holder inequality to (4.4) we find 



(u e ® u e - v <8> v + I)/ e dv 



dx ds (4.22) 



<[([ I r\v\ 2 + f> e \ 2 dvdx) 2 D^nUs < c ( [ D^n 

Jo \Jm d JR d / \Jo 

where the last inequality follows from the entropy bound (Proposition 2.1) and 

— 



e |«T = / fvu e dv<[ I /V dvj (J f\u e \ z dv 

1 

We conclude by combining (4.21) and (4.22). □ 
4.5. (4) The last term. Finally, we have: 

Lemma 4.9. Assuming that U = (g, gu) and V — (q,qv) are such that u € 
L° a (M. d ), q, g E L 1 (M. d ), there exists a constant C such that 

K(x,y)q(x)(g(y) - q(y))[u(y) - u(x)][v(x) ~ u(x)] dxdy 

, (4-23) 
<C||u||L~(||g|Ui+||g||ii) / £ ( v \ u ) dx 

JTS. d 

Proof. We have 

K(x, y)q(x)(g(y) - q(y))[u(y) - u(x)][v(x) - u(x)} dxdy 

1 fl 

m&x{q(y),g(y)} \v[x) - u(x)\ dxdy 

- 2Ml ~ {L L K{Xl v)q[x) min { W) ' W) } {e{v) q{v)f dxdy ) 

■ 1/2 

K(x, y)q{x) max {q(y), g(y)} \v(x) - u(x)\ 2 dxdy 
And so using (4.3) and the fact that K(x,y) < C, we deduce: 

K(x, y)q(x)(g(y) - q(y))[u(y) - u(x)][v(x) - u(x)} dxdy 

< C||u|| z -||g|lil a (||g|Ui + \\q\W) 1 ' 2 Up(q\Q)(y) dy^'* 

x (yj q(x)[v(x) - u(x)] 2 dxj 

< C||u|U-||g||£/ 2 (|M| £ i + |MU0 1/2 / £(Y\U) dx. 

JR d 

a 
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4.6. Proof of Proposition 4.1. We recall that 

£(£/ e (0)|J7(0)) = 0. 

For any t £ (0,t), we integrate (4.18) over (0,i) and apply Lemmas 4.7, 4.8 and 
4.9, and , to obtain the inequality 

£{U c \U){t) dx 

>. 

t rt r f 

K(x,y)g e (x)g e (y) [(u e (x) - u(x)) - (u e (y) - u(y))} 2 dx dy ds 
f f d t E{W) + P e V x $ + \ f K(x 1 y)g^x)g i (y)[u t (x)-u t (y)} 2 dydxds 



la Jm d 



< 



+ CV^ + C f f £(U e \U)dxdt. 

JO JR d 

Lemma 4.6 now implies: 

f £{U e \U)(t) dx 
Jm d 

+ 111 I K(x, y)g e (x)g<(y) [(u*(x) - u(x)) - {u\y) - u{y))] 2 dxdydt 
A Jo Jm d Jm d 

<Ce I T(f e (s))ds + C'V^ + C f f £(U e \U) dxdt. (4.24) 
Jo Jo Jwi d 



and (2.10) implies 

We deduce 

f £(U e \U)(t) dx 

Jwi d 

I/"*/"/ 



F(f{s))ds<e 



Cat 



\ \ K{x, y)g e (x) Q e {y) [(u e (x) - u(x)) - (u e (y) - u(y))f dxdydt 

- -jki JS. d JM d 

<C(T)Ve + c/ f £{U e \U) dxdt. (4.25) 

JO J«. d 

which completes the proof of Proposition 4.1. 
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Appendix A. Local well-posedness of the Euler-flocking system 

The purpose of this appendix is to prove the existence of a local-in-time unique 
smooth solution of the Euler-flocking equations. In particular, the objective is to 
prove Theorem 2.4 which we relied upon to prove our main result (Theorem 3.1). 

First, we observe that the system (1.5)-(1.6) is a 4 x 4 system of conservation 
laws which can be written in the following equivalent form when the solution is 
smooth: 

\d t g + uS7 x g + gdiv x u = 

\g(d t u + u\7 x u) +S7 x g = F(g,u,\7 x <i>). 

Here, 

F(q,u,V$(x))= [ K{x,y)g(x)g{y)[u(y)-u{x)} dy ~ g\7 x <S>. (A.2) 

Js. d 

In what follows, we shall need the Sobolev space H s (M. d ) given by the norm 

hWl = £ / \D a g\ 2 dx. 

\ol\<s 

For g e L°°([0,T];H s ), define 

INIkr = sup ||ff(-,t)|| s . 

0<t<T 

Now, consider the Cauchy problem of (A.l) with smooth initial data: 

(£>,tt)| t= o = {g ,u )(x). (A.3) 

The values of the vector 

V 



u, 

lie in the state space G, which is an open set in K 4 . The state space Q is introduced 
because certain physical quantities such as density should be positive. Indeed, by 
invoking the method of characteristics for the continuity equation the following 
lemma holds: 

Lemma A.l. If (g,u) G C 1 (M 3 x [0,T]) is a uniformly bounded solution of (1.5) 
with g{x,0) > 0, then g(x,t) > on R 3 x [0,t\. 



System (A.l) can be written in the form 

dt ( Q ) + uV x 
or equivalently 



g\ ( g\ ( g divj; u 



Q\ V7 ( P \ ( 



dt y +V X - [ES£ +i oggh ) - {l F (g,u,WMx)) ( A - 4 ) 

with P = gu and F given in (A.2). 

It turns out that (A.4) has the following structure of symmetric hyperbolic sys- 
tems: For all w € Q, there is a positive definite matrix Aq(w) that is smooth in w 
and satisfies: 

Co % < A {w) < coh, 
with a constant Cq uniform for w G Q\ C Q\ C Q such that 

Ai(w) = A (w)Vfi(w) 

is symmetric. 

Here, Vfi{w),i = 1,2,3 are the 4x4 Jacobian matrices and I4 is the 4x4 
identity matrix. 
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The matrix 

is called the symmetrizing matrix of system (A. 4). Multiplying (A. 4) by Aq we 
obtain 

A (w)d t w + A(w)V x w = G(w, V x $) (A.5) 
with smooth initial data 

w o =w(x,0). (A.6) 
Here, A(w) = (Ai(w), . . . , As(w)) denotes a matrix whose columns 
Ai{w) — Ao(w)V x /i(w) are symmetric, and 

P 



fH 



tr 

with P = qu and V x fi(w) : are the 4x4 Jacobian matrices and I3 denotes the 3x3 
identity matrix. The 4x1 vector G in (A.5) is given 

G ( W , v.*) = M«) (i F{Q , ° v,«)) - U, u° v.*)) (AJ) 

which implies that 

We are now ready to prove the local existence of smooth solutions. 

Theorem A. 2. Assume w = (go, u ) G H s n L°°(IR 3 ) with s > 5/2 and qq(x) > 
and that V^4> 6 _ff s . TTien there is a finite time T G (0, 00), depending on the H s 
and L°° norms of the initial data, such that the Cauchy problem (1.5)-(1.6) and 
(A. 3) has a unique bounded smooth solution w = (g, u) G C 1 (M 3 x [0, T]), with g > 
for all (x,t) £l 3 x [0,T], and (g,u) G C([0,t];H s ) n C x ([0, T]; if 5 " 1 ). 

Theorem A. 2 is a consequence of the following theorem on the local existence of 
smooth solutions, with the specific state space Q — {(g, u) T : g > 0} C R 4 for the 
inhomogeneous system (A.5). 

Theorem A. 3. Assume that wq : K d — > Q is in H s n L°° with s > § + 1. Then, for 
the Cauchy problem (A.5)-(A.6) ; there exists a finite timeT = T(\\wq\\ s ,\\wo\\l°°) G 
(0,oo) such that there is a classical solution w G C 1 (M 3 x [0,T]) with w(x,t) G Q 
for(x,t) eR d x [0,T] and we C([0, T]; iJ s ) n C 1 ([0, T];H S ^ 1 ). 

Proof. The proof of this theorem proceeds via a classical iteration scheme method. 
An outline of the proof of Theorem A. 3 (and therefore Theorem A. 2) is given below. 
Consider the standard mollifier 

n(x)eC^(R 3 ), supp rj(x) C {x; \x\ < 1}, T](x)>0, / J7(as)dx = 1, 
and set 

r\e = e~ d n(x/e). 
Define the initial data w% G C°°(]R 3 ) by 



™a( x ) = Ve h *w (x) = / r) ek (x-y)w (y)dy, 

JM. 3 

where = 2~ k eo with eo > constant. We construct the solution of (A.5)-(A.6) 
using the following iteration scheme. Set w°(x,t) = Wq(x) and define w k+1 (x,t), 
for k = 0, 2, . . . , inductively as the solutions of linear equations: 

j A {w k )d t w k+1 + A(w k )V x w k+1 = G(w k , V*$), 

K +1 | t =o = ^ +1 (z)- 
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By well-known properties of mollifiers it is clear that: 

\\wq - w \\ s 0, asfc -> oo, and \\w$ - w \\ < C e fe ||w ||i, 

for some constant Cq. Moreover, w k+1 £ C ca (M. d x [0,Tfc]) is well-defined on the 
time interval [0,71-], where T k > denotes the largest time for which the estimate 
\\\w h — Wo||| Sj T fc < C\ holds. We can then assert the existence of T„ > such that 
Tk > T* (To = oo) for k = 0, 1, 2, . . . , from the following estimates: 

lh fc+1 - <\\\s,t, < C x , \H +1 \\ S ^ T , < C 2 , (A.9) 

for all k = 0, 1, 2, ... , for some constant C2 > 0. From (A. 8), we get 

A Q (w k )d t (w k+1 - w k ) + A(w k )V(w k+1 - w k ) = E k + G k , (A.10) 

where 

f E k = -(A (w k )-A (w k - 1 )d t w k -(A(w k )-A(w k - 1 ))Vw k ) 
\ G k = G(w k ,V x $). 
From the standard energy estimate for the linearized problem (A. 10) we get 

\\\w k+1 - w k \\\ ,T < Ce CT (\\w k+1 - w k \\ + T\lE k \\ ,T + T||G fc |l 0l T). 
Taking into consideration the property of mollification, relation (A.9), we have that 
K +1 - ug||o < C2- k , \\E k \\ 0tT < C\\\w k - w k -% iT . 

Note that, 

||G fc |lo,T= sup ||G fc (.,t)IU 

0<t<T 



sup 

0<t<T 



K(x,y) Q k (x)g k (y)[u k (y) - u k (x)] dy - g k V x $(x) 
Vu k (y)\ 2 



Hx)\ 2 



<M[ \vVv k \ 2 dy + C = M [ |VM fc . ( f dy + C. 
Here we use the well known result that states: For any v £ J7 1 (M 3 

r dx < M \Vv\ 2 dx. 



1 

For small T such that C 2 T exp CT < 1 one obtains 

00 

£|||™ fc+1 -™ fc |||o,r<oo, 

k=l 

which implies that there exists w £ C([0, T]; L 2 (M. 3 )) such that 

lim |||u) fc - w\\\ 0T = 0. (A.12) 

k— >oo 

From (A.9), we have |||w fe ||| s ,t + \\\w k ||| s -i,t < C, and w k (x,t) belongs to a bounded 
set of G for (x, t) £ K d x [0, T]. Then by interpolation we have that for any r with 
< r < s, 

\\\w k - w%, T < C s \\\w k - v l \\\l^ /a \\\w k - w l \\\:^ < \\\w k ~ w%^ r/s . (A.13) 
From (A.12) and (A.13), 

,k 

— t 

for any < r < s. Therefore, choosing r > I + 1, Sobolev's lemma implies 

w k -> w in C([0,i];C ,1 (M 3 )). (A.14) 

From (A.10) and (A.14) one can conclude that w k -> w in C([0, T]; C(R 3 ), w £ 
C 1 (R 3 x [0, T]), and w(x,t) is a smooth solution of (A.5)-(A.6). To prove u £ 



lim \\\w - w\\ r .,T = 

k— >oo 
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C([0,T];H s ) n ^([OjT];^" 1 ), it is sufficient to prove u £ C{[0,T];H s ) since it 
follows from the equations in (A.5) that u £ C^fO, T];H s ~ l ). □ 

Remark A. 4. The proof of Theorem A. 2 follows the line of argument presented 
by Majda [22], which relies solely on the elementary linear existence theory for 
symmetric hyperbolic systems with smooth coefficients (see also Courant-Hilbert 

[8])- 

Lemma A.5. Let (g, u) be a sufficiently smooth solutions to (1.5)-(1.6). The energy 
E(u) satisfy the following entropy equality 

dt I E{u)dx+\f K(x,y)g(x)g(y)\u(x)-u(y)\ 2 dydx = 0. (A.15) 

JR d 1 JR d 

Proof. The result is obtained by the following standard process: first we multiply 
the continuity equation (1.7) by (glogg)' and the momentum equation (1.8) by 
the velocity field u, and we add the resulting relations. Next we integrate over the 
domain taking into account that (g,u) is sufficiently smooth and using that 

d t / g<&dx= / d t g&dx = — / div(gu)§ dx = / / guV x <&dx. 

JR d JR d JR d JR d JR d 

/ u(x) K(x,y)g(x)g(y)[u(y) - u(x)\dydx 

JM d JR d 



and 



K(x,y)g(x)g(y)\u(x) -u(y)\ dydx 



2 

we arrive at (A.15). The result for a solution (g,u) with the regularity established 
in Theorem A. 2 is established using a density argument. □ 



A.l. Proof of Theorem 2.4. The proof of Theorem A. 2 can be further reduced to 
verifying that u(x, t) is strongly right continuous at t = 0, since the same argument 
works for the strong right-continuity at any other t £ (0, T] and the strong right- 
continuity on [0, T) implies the strong left-continuity on (0, T] because the equations 
in A.5 are reversible in time. Taking this into consideration estimate (A.15) implies 
the result. □ 
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